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Abstract: We give an interpretation of the t = 1 specialization of the 
modified Macdonald polynomial as a generating function of the energy 
statistics defined on the set of paths arising in the context of Box-Ball 
Systems (BBS-paths for short). We also introduce one parameter gen- 
eralizations of the energy statistics on the set of BBS-paths which all, 
conjecturally, have the same distribution. 

Resume: Nous donnons une interpretation de la specialisation a t = 1 
du polynome de Macdonald modifie comme fonction generatrice des 
statistiques d'energie definies sur I'ensemble des chemins qui apparais- 
sent dans la theorie des Systemes BBS (BBS-chemins). Nous presentons 
egalement des generalisations a un parametre de la statistique d'energie 
sur les chemins BBS qui toutes, conjecturalement, ont la meme distribu- 
tion. 
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1 Introduction 



The purpose of the present paper is two-fold. First of all we would like to draw atten- 
tion to a rich combinatorics hidden behind the dynamics of Box-Ball Systems, and 
secondly, to connect the former with the theory of modified Macdonald polynomials. 
More specifically, our final goal is to give an interpretation of the Kostka-Macdonald 
polynomials Kx^^{q,t) as a refined partition function of a certain box-ball systems 
depending on initial data A and /i. 

Box-Ball Systems (BBS for short) were invented by Takahashi-Satsuma [29| [28] 
as a wide class of discrete integrable soliton systems. In the simplest case, BBS 
are described by simple combinatorial procedures using boxes and balls. One can 
see the simplest but still very interesting examples of the BBS by the free software 
available at [26]. Despite its simple outlook, it is known that the BBS have various 
remarkably deep properties: 

• Local time evolution rule of the BBS coincides with the isomorphism of the 
crystal bases [71[2]. Thus the BBS possesses quantum integrability. 

• BBS are ultradiscrete (or tropical) limit of the usual soliton systems [30| 120]. 
Thus the BBS possesses classical integrability at the same time. 

• Inverse scattering formalism of the BBS [19] coincides with the rigged con- 
figuration bijection originating in completeness problem of the Bethe states 
[UllIE], see also [25]. 

Let us say a few words about the main results of this note. 

• We will identify the space of states of a BBS with the corresponding weight 
subspace in the tensor product of fundamental (or rectangular) representations 
of the Lie algebra Qi{n). 

• In the case of statistics tau, our main result can be formulated as a computation 
of the corresponding partition function for the BBS in terms of the values of 
the Kostka-Macdonald polynomials at t = 1. 

• In the case of the statistics energy, our result can be formulated as an inter- 
pretation of the corresponding partition function for the BBS as the g-weight 
multiplicity of a certain irreducible representation of the Lie algebra Q\{n) in 
the tensor product of the fundamental representations. We expect that the 
same statement is valid for the BBS corresponding to the tensor product of 
rectangular representations. 

Let us remind that a q- analogue of the multiplicity of a highest weight A in 
the tensor product (S)a=i ^sa-^ra highest weight Sa ojr^, a = 1,...,L, 

irreducible representations Kaw^o of the Lie algebra fll(n) is defined as 

L 

g-Mult [V, : (g) K„^,J = J2 K,,R K,,x{q), 

a=l rj 
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where K^^^r stands for the parabohc Kostka number corresponding to the se- 
quence of rectangles R := {(s^")}a=i,...,L, see e.g. [15], [T8] . 

A combinatorial description of the modified Macdonald polynomials has been 
obtained by Haglund-Haiman-Loehr [5] . In Section [5] we give an interpretation of 
two Haglund's statistics in the context of the box-ball systems, i.e., in terms of the 
BBS-paths. Namely, we identify the set of BBS paths of weight a with the set V{a) 
which is the weight a component in the tensor product of crystals corresponding to 
vector representations. We have observed that from the proof given in [5] one can 
prove the following identity 

J2 ginv,(p)^maj,(p) ^ ^ K^ J{^,^(q, t) , (l) 

see Proposition 16.21 and Corollary 16. 3[ One of the main problems we are interested 
in is to generalize the identity Eq.([T]) on more wider set of the BBS-paths. 

Our result about connections of the energy partition functions for BBS and q- 
weight multiplicities suggests a deep hidden connections between partition functions 
for the BBS and characters of the Demazure modules, solutions to the g-difference 
Toda equations, cf.[3], ... . 

As an interesting open problem we want to give raise a question about an in- 
terpretation of the sums -^t?,_r -^»;,a(?,^), where Kr,^\{q,t) denotes the Kostka- 
Macdonald polynomials [21], as refined partition functions for the BBS correspond- 
ing to the tensor product of rectangular representations R = {(s^")}i<a<n- In other 
words, one can ask: what is a meaning of the second statistics (see [5]) in the 
Kashiwara theory [11] of crystal bases (of type A) ? 

This paper is abbreviated and updated version of our paper [17]. The main 
novelty of the present paper is the definition of a one parameter family of statistics 
on the set of BBS-paths which generalizes those introduced in [17], see Conjecture 
17.21 It conjecturally gives a new family of MacMahonian statistics on the set of 
transportation matrices, see [T5] . 

Organization of the present paper is as follows. In Section [2] we outlook the basic 
definitions and facts related to the Kashiwara's theory of crystal base in the case of 
type An\ We also remind definitions of the combinatorial i?-matrix and definition 
of the energy function. We illustrate definitions by simple example. In Section [3l 
we introduce the energy statistics and the set of the BBS. In Section H] we remind 
definition of box-ball systems and state some of their simplest properties. In Section 
O we remind definition of the Haglund's statistics and give their interpretation in 
terms of the BBS-paths. Sections [6] and [7] contain our main results and conjectures. 
In particular it is not difficult to see that Haglund's statistics maj^ and inv^ do not 
compatible with the Kostka-Macdonald polynomials for general partitions A and /i. 
In Section Owe state a conjecture which describes the all pairs of partitions (A,/i) 
for those the restriction of the Haglund-Haiman-Loehr formula on the set of highest 
weight paths of shape /i coincide with the Kostka-Macdonald polynomial Kx^^{q, t). 



3 



2 Kirillov— Reshetikhin crystal 

2.1 An'' type crystal 

Let wt^ be a U'g{g) Kirillov-Reshetikhin module, where we shall consider the case 
g = An^ . The module Ws^^ is indexed by a Dynkin node r E I = {l,2,...,n} 
and s G Z>o- As a -module, wi^"* is isomorphic to the irreducible module 

corresponding to the partition (s^). For arbitrary r and s, the module Wt^ is known 
to have crystal bases [HI [10], which we denote by B""'^. As the set, B^'^ is consisting 
of all column strict semi-standard Young tableaux of depth r and width s over the 
alphabet {1, 2, . . . , r?, + 1}. 

For the algebra An, let P be the weight lattice, {Aj G P\i G /} be the funda- 
mental roots, {oj G P\i G /} be the simple roots, and {hi G B.omz{P,'Z)\i G /} be 
the simple coroots. As a type An crystal, B = B^'^ is equipped with the Kashiwara 
operators e^, : B — y B U {0} and wt : B — y P (i E I) satisfying 

f,{b) = b' ^ i,{b') = b ifb,b'eB, 

wt{m) = wt{b)-ai am^B, 

(/i„wt(6)) = ipi{b) -e^{b). 

Here (■, ■) is the natural pairing and we set ei{b) = max{m > | e™6 ^ 0} and 
(pi{b) = max{m > | f^b ^ 0}. Actions of the Kashiwara operators Cj, fi for z G / 
coincide with the one described in [12] . Since we do not use explicit forms of these 
operators, we omit the details. See [23] for complements of this section. Note that 
in our case An, we have P = Z""*"^ and = — e^+i where is the i-th canonical 
unit vector of We also remark that wt(b) = (Ai, ■ ■ ■ , A^+i) is the weight of b, 

i.e., Aj counts the number of letters i contained in tableau b. 

For two crystals B and B', one can define the tensor product B B' = {b ® b' \ 
b E B,b' E B'}. The actions of the Kashiwara operators on tensor product have 
simple form. Namely, the operators Cj, fi act o\i B ® B' by 



ei{b®b') 
kb®b') 



eib®b' if ipi{b) > 6i{b') 

b®eib' li ipi{b) <ei{b'), 

~fib®b' ii ^iib) > Eiib') 

b®~fib' if (^,(6) < 6,(6'), 



and wt(6 ® b') = wt(6) + wt{b'). We assume that (8> 6' and 6 (8> as 0. Then it is 
known that there is the unique crystal isomorphism R : B'^'^ (g> B'^''^' ^ B''''^' ® B'^'^. 
We call this map (classical) combinatorial R and usually write the map R simply 
by ~. 

Let us consider the affinization of the crystal B. As the set, it is 

Aff (5) = {b[d] \beB,deZ}. (2) 

There is also explicit algorithm for actions of the affine Kashiwara operators eo, 
/o in terms of the promotion operator [27]. For the tensor product b[d\ ® b'[d'] G 
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AfF(i?) (8) Aff (i?'), we can lift the (classical) combinatorial R to affine CRSG clS follows: 
h[d] ® h'[d'] ~ h'[d' - H{h ® h')] ® b[d + H{b ® b')], (3) 

where 6 ® 6' ~ 6' ® 6 is the isomorphism of (classical) combinatorial R. The function 
H{b®b') is called the energy function and defined by a certain set of axioms. We will 
give explicit forms of the combinatorial R and energy function in the next section. 

2.2 Combinatorial R and energy function 

We give an explicit description of the combinatorial i?-matrix (combinatorial R for 
short) and energy function on B''''^®B^' . To begin with we define few terminologies 
about Young tableaux. Denote rows of a Young tableaux Y by yi, 1/2, ■■ - yr from top 
to bottom. Then row word row(Y) is defined by concatenating rows as row{Y) = 
yryr~i ■ ■ - Hi- Let X = {xi, X2, ■ ■ .) and y = {yi, 2/2, . . .) be two partitions. We define 
concatenation of x and y by the partition {xi + yi, X2 + y2, ■ ■ ■)■ 

Proposition 2.1 (t27j) b b' e B''^' ® 5^'-"' is mapped tob'(^be B''^'' ® 5'^'" 
under the combinatorial R, i.e., 

6 (g) fe' ~ 6' (g) 6, (4) 

ij and only if 

{b' ^row{b)) = (b^row(b')). (5) 

Moreover, the energy function H{b ® b') is given by the number of nodes of {b' ^ 
row{b)) outside the concatenation of partitions {.$''') and {s'^ ). 

For special cases of B^'^ (g) B^'^' , the function H is called unwinding number in 
[22] . Explicit values for the case b^b' G B^'^^B^'^ are given by H{b®b') = x{b < b') 
where x(True) = 1 and x(False) = 0. 

In order to describe the algorithm for finding b and b' from the data {b' <— row{b)), 
we introduce a terminology. Let y be a tableau, and Y' be a subset of Y such that 
Y' is also a tableau. Consider the set theoretic subtraction 9 = Y\Y'. If the number 
of nodes contained in 61 is r and if the number of nodes of 9 contained in each row 
is always or 1, then 9 is called vertical r-strip. 

Given a tableau Y = {b' ^ row{b)), let Y' be the upper left part of Y whose 
shape is (s^). We assign numbers from 1 to r's' for each node contained in ^ = F \y 
by the following procedure. Let 6*1 be the vertical r'-strip of 9 as upper as possible. 
For each node in ^1, we assign numbers 1 through r' from the bottom to top. Next 
we consider 9\9i, and find the vertical r' strip ^2 by the same way. Continue this 
procedure until all nodes of 9 are assigned numbers up to r's'. Then we apply inverse 
bumping procedure according to the labeling of nodes in 9. Denote by ui the integer 
which is ejected when we apply inverse bumping procedure starting from the node 
with label 1. Denote by Yi the tableau such that (Yi ^ ui) = Y. Next we apply 
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inverse bumping procedure starting from the node of Yi labeled by 2, and obtain 
the integer U2 and tableau Y2. We do this procedure until we obtain Ur's' and Yr's'- 
Finally, we have 

b' = ^ Ur's'Ur's'~l ■ ■ ■ Ui), b = Yr's'- (6) 



Example 2.2 Consider the following tensor product: 



1 


1 


4 


2 


3 


6 



2 


3 


3 


4 


4 


5 



From fe, we have row{b) = 236114, hence we have 





2 


3 




1 


3 


4 


^ 236114 




4 


5 





1 


1 


3 


43 


2 


2 


6 




36 


32 




45 


4i 




54 







Here subscripts of each node indicate the order of inverse bumping procedure. For 
example, we start from the node 4i and obtain 





1 


2 


3 


4 


\ 


1 


1 


3 


4 


1 


2 


3 


43 




2 


3 


6 




2 


2 


6 




2 


3 


6 






3 


4 




^ 1 = 


= 3 


3 




, therefore, Yi = 


36 


42 






4 








4 


4 






45 






v 


5 






/ 


5 








54 







Ml = 1. 



Next we start from the node 42 of Yi. Continuing in this way, we obtain ^6^5 ■ ■ - Ui 

Since (0 ^ 321421) = | 2 I 2 I, we obtain 



321421 and K; 



3 


3 


4 


4 


5 


6 



1 


1 


4 


2 


3 


6 



2 


3 




1 


1 


3 


4 




2 


2 


4 


5 




3 


4 



3 


3 


4 


4 


5 


6 



H 



1 


1 


4 


2 


3 


6 



Note that the energy function is derived from the concatenation of shapes of b and 
6', i.e. 
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3 Energy statistics and its generalizations on the 
set of paths 

For a path 61 (g) ^2 ® ■ ■ ■ ® &l e -B'"'''" ® -B''''*'' ® ■ ■ ■ ® S^'-f-'^-^', let us define elements 
fej*"* G for z < J by the following isomorphisms of the combinatorial R; 

61 ® ^2 ® 
~ 61 (g) 62 ® 

^ 61 ® 62 ® 

where we have written bk ® l)-^^^'' ~ ® b'f^ assuming that bj'^ = bj. 
Define the statistics maj(p) by 

mai{p) = J2H{k^bf+'^). (8) 

i<j 

For example, consider a path a = ai (g) 02 ® ■ ■ ■ ® G {B^'^)'^^. In this case, we 
have ttj"^ = tti, since the combinatorial R act on B^'^ (g) B^'^ as identity. Therefore, 
we have 

L-l 

maj(a) = J^(L - < fli+i)- (9) 

Define another statistics tau as follows. 
Definition 3.1 For the path p E B""^''^ B''^''^ O ■ ■ ■ (g) S^'-^-^^, c?e/ine r^'" 

r'''''(p) = maj(w^'') (gip), (10) 
where Us is the highest element of B'''''^. 

Here the highest element ui^^ G -B*^'* is the tableau whose i-th row is occupied by 

. In particular, the statistics r*"'^ on B^'^ 
type paths a G (^B^'^)^^ has the following form; 

L-l 

T''\a) = L-x{r< ai) + - i)x{ai < Oi+i), (11) 

i=l 

where ai denotes the first letter of the path a. Note that r^'^ is a special case of the 
tau functions for the box-ball systems [201 [21] which originates as an ultradiscrete 
limit of the tau functions for the KP hierarchy [9] . 

Definition 3.2 For composition fi = (/ii,/i2,--- ,/in), write fi[ij = J2]=if^j ^^^^ 
convention ^[0] = 0. Then we define a generalization of t^'^ by 

n 

r;'\a) = J2r'-'\a^^), (12) 

1=1 
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(g) bi 
bi-i ® bi 



6' 



(7) 



integers i. For example, = 



1 


1 


1 


1 


2 


2 


2 


2 


3 


3 


3 


3 



where 

«w = ® «M[,-i,+2 ® ■ ■ ■ ® e (5^'^)®'^^ . (13) 

Note that we have a = a[i] a[2] ■ ■ ■ ® a[n], i-e., the path a is partitioned according 
to /i. 



4 Box-ball system 

In this section, we summarize basic facts about the box-ball system in order to 
explain physical origin of r^'^. For our purpose, it is convenient to express the 
isomorphism of the combinatorial R: a® 6 ~ 6'(x3a' by the following vertex diagram: 



a- 



y 

Successive applications of the combinatorial R is depicted by concatenating these 
vertices. 



Following [71 [2] , we define time evolution of the box-ball system t/"'' . Let 



G -B"'' be the highest element and hi G B^''^\ Define uj"'' and b'^ G B^""^^ by the 
following diagram. 



u 



(a) 
1,0 



(a) 



u 



(a) 
LO' 



■U 



(a) 
LI 



■u 



(a) 
1,2 ' 



U 



(a) 

l,L-l' 



b'o 



■ U 



(a) 
l,L 



b'r 



(14) 



,(«) ._ 



u["'j are usually called carrier and we set u]"q 
by 

T[''\b) = b' = b[^b'^ 



U 



(a) 



Then we define operator T, 



(a) 



h'L- 



•2yy ■ ■ ■ '-^ ^L- (15) 

Recently [25], operators T^""' have used to derive crystal theoretical meaning of the 
rigged configuration bijection. 

It is known ([19] Theorem 2.7) that there exists some / G Zi>o such that 



T, 



(a) 



rp{a) _ rji(a) 



(16) 



If the corresponding path is 6 G (5^'^)®^, we have the following combinatorial 
description of the box-ball system [291 128]. We regard [T] G B^'^ as an empty box of 
capacity 1, and |T| G B^'^ as a ball of label (or internal degree of freedom) i contained 
in the box. Then we have: 

Proposition 4.1 ([?]) For a path b G (S^'^)®-^ of type A^n \ T^\b) is given by the 
following procedure. 



1. Move every hall only once. 
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2. Move the leftmost ball with label n + 1 to the nearest right empty box. 



3. Move the leftmost ball with label n+1 among the rest to its nearest right empty 
box. 

4. Repeat this procedure until all of the balls with label n + 1 are moved. 

5. Do the same procedure 2-4 for the balls with label n. 

6. Repeat this procedure successively until all of the balls with label 2 are moved. 

There are extensions of this box and ball algorithm corresponding to generalizations 
of the box-ball systems with respect to each affine Lie algebra, see e.g., [8]. Using this 
box and ball interpretation, our statistics T^'^(fe) admits the following interpretation. 



Theorem 4.2 ([20] Theorem 7.4) For a path h e (fii-i)®-^ of type A 



(1) 



r 



1,1 



coincides with number of all balls 2, ■ ■ ■ ,n + 1 contained in paths b, T^\h), 



(1)\L-1 



(b). 



Example 4.3 Consider the path p = a (g) 6 where a = 4311211111, b = 4321111111. 
Note that we omit all frames of tableaux of B^'^ and symbols for tensor product. We 
compute T(io,io)(p) by using Theorem 14.21 According to Proposition 14.11 the time 
evolutions of the paths a and b follows: 



4 


3 


1 


1 


2 


1 


1 


1 


1 


1 




4 


3 


2 


1 


1 


1 


1 


1 


1 


1 


1 


1 


4 


3 


1 


2 


1 


1 


1 


1 




1 


1 


1 


4 


3 


2 


1 


1 


1 


1 


1 


1 


1 


1 


4 


1 


3 


2 


1 


1 




1 


1 


1 


1 


1 


1 


4 


3 


2 


1 


1 


1 


1 


1 


1 


4 


1 


1 


3 


2 




1 


1 


1 


1 


1 


1 


1 


1 


1 


4 


1 


1 


1 


1 


1 


1 


4 


1 


1 


1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


4 


1 


1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


4 


1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


4 




1 


1 


1 


1 


1 


1 


1 


1 


1 


1 



Here the left and right tables correspond to a and b, respectively. Rows of left (resp. 
right) table represent a, T^\a), {T^^)^{a) (resp., those for b) from top to 
bottom. Counting letters 2, 3 and 4 in each table, we have T^'\a) = 16, T^'^ib) = 10 
and we get T(\qiq)(p) = 16 + 10 = 26, which coincides with the computation by 
Eg. (ITT]) . Meanings of the above two dynamics corresponding to paths a and b are 
summarized as follows: 

(a) Dynamics of the path a. In the first two rows, there are two solitons (length 
two soliton 43 and length one soliton 2), and in the lower rows, there are also 
two solitons (length one soliton 4 and length two soliton 32). This is scattering 
of two solitons. After the scattering, soliton 4 propagates at velocity one and 
soliton 32 propagates at velocity two without scattering. 

(6) Dynamics of the path b. This shows free propagation of one soliton of length 
three 432 at velocity three. 
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5 Haglund's statistics 



Tableaux language description For a given path a = ai ® a2 ® ■ ■ ■ ® ai G 

(5^'^)®^, associate tabloid t of shape /i whose reading word coincides with a. For 
example, to path p = abcdefgh and the composition jj, = (3,2,3) one associates the 
tabloid 



c 


b 


a 


e 


d 




h 


9 


f 



(17) 



Denote the cell at the i-th row, j-th column (we denote the coordinate by 
of the tabloid t by tij. Attacking region of the cell at is all cells {i,k) with 
k < j OT {i + l,k) with k > j. In the following diagram, gray zonal regions are the 
attacking regions of the cell 




Follow [5], define |Invjj| by 

|lnvjj| = j^{{k,l) G attacking region for \ tki > tij}. (18) 
Then we define 

|lnv^(a)| = J2 I^^^^J-I- (19) 
If we have < tij, then the cell is called by descent. Then define 

Des^(a)= Yl if'^-j)- (20) 

all descent 

Note that (/ij — j) is the arm length of the cell 

Path language description Consider two paths a(i),a(2) e (^I'l)®'^. We denote 
by a^^^ ® a^^-* = ai ® 02 ® ■ ■ ■ ® a2^. Then we define 

Inv(^,^)(a(i\a(2)) = J2Y1 ^("'^ < (21) 

k=l i=k+l 

For more general cases a^^^ E ^^d a*^^^ G (iJi'i)®'^2 satisfying fii > fj,2, we 

define 

lnv(,,,,,)(a«,a(2)) := Inv(^,^,)(a«, 1®(^-^^) ® a^^)). (22) 
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Then the above definition of |Inv^(a)| is equivalent to 

n-l 

|Inv^(a)| = ^Inv(^^,^^^^). (23) 

Consider two paths a^^^ G {B^'^)^^'^ and a^^^ G {B^'^)®^'^ satisfying /xi > //a- 
Denote a = a^^^ ® a*-^''. Then define 

Ml 

Des(^i,/,2)(«) = XI - (^1 - ^2) - l)x(afc < Ofc+^J. (24) 

fc=Mi-M2+l 

For the tableau T of shape corresponding to the path a, we define 

n 

Des/,(r) = J^Des(/,^,^^^,)(a[i] ® a^^+i]). (25) 
Definition 5.1 ([4]) For a path a, statistics maj^ zs defined by 

Ml 

maj^(a) = Xmaj(fi,i (g) ^2,1 ® ■ ■ ■ ® t;,^,*). (26) 
1=1 

and inv^(a) is defined by 

inv^(a) = |Inv^(a)| - Des^(a). (27) 

If we associate to a given path p G V{X) with the shape // tabloid T, we sometimes 
write mai^{p) = maj(T) and inv^(p) = inv(T). 

6 Haglund— Haiman— Loehr formula 

Let Hfj^{x;q,t) be the (integral form) modified Macdonald polynomials where x 
stands for infinitely many variables Xi,X2, - ■ ■ ■ Here H^{x; q, t) is obtained by simple 
plethystic substitution (see, e.g., section 2 of |6]) from the original definition of the 
Macdonald polynomials [21]. Schur function expansion of H^{x] q,t) is given by 

H^{x;q,t) = J2KxA(l^t)sx{x), (28) 

A 

where Kx^q^t) stands for the following transformation of the Kostka-Macdonald 
polynomials: 

Here we have used notation n{ij) = ~ l)/^*- Then the celebrated Haglund- 

Haiman-Loehr (HHL) formula is as follows. 
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Theorem 6.1 ([5]) Let cr : /i — )■ Z>o be the filling of the Young diagram fi by 
positive integers Z>o, and define = Yluefi^^^iu)- Then the Macdonald polynomial 
Hfj^{x;q,t) have the following explicit formula: 

H^{x-q,t)= g"^"(")r^j('^)x". (30) 

From the HHL formula, we can show the foUowing formula. 
Proposition 6.2 For any partition yU and composition a of the same size, one has 

E r-'^^'^'t^^'-^'^ = K,,aK,A<l, t), (31) 

where V{a) stands for the set of type B^'^ paths of weight a = (ai,a2, • • • ,an+i) 
and 7] runs over all partitions of size \^\. 

Corollary 6.3 The (modified) Macdonald polynomial H^{x\q,t) have the following 
expansion in terms of the monomial symmetric functions m\[x): 

H,{x-q,t) = Y.\Y. ^a(x), (32) 

\pGP(A) / 

where A runs over all partitions of size |/x|. 

To find combinatorial interpretation of the Kostka-Macdonald polynomials Kx^^{q, t) 
remains significant open problem. Among many important partial results about 
this problem, we would like to mention the following theorem also due to Haglund- 
Haiman-Loehr: 

Theorem 6.4 (|5J Proposition 9.2) // /ii < 2, we have 

Kx,f,iq,t)= J2 ^'""''(P^r^j^^P), (33) 

where V+IX) is the set of all highest weight elements ofV{X) according to the reading 
order explained in Eg. [FID . 

It is interesting to compare this formula with the formula obtained by S. Fishel [1], 
see also [H], [18]. 

Concerning validity of the formula Eq.( l33|) . we state the following conjecture. 

Conjecture 6.5 Explicit formula for the Kostka-Macdonald polynomials 

Kx,^,{q,t)= (34) 

is valid if and only if at least one of the following two conditions is satisfied, 
(i) /t^i ^ 3 and < 2. 
(ii) X is a hook shape. 
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7 Generating function of tau functions 



In [T7], we give an elementary proof for special case t = 1 of the formula Eq.f l3ip in 
the following form. 

Theorem 7.1 Let a be a composition and be a partition of the same size. Then, 



^ K^^^,{q, 1). 



Conjecture 7.2 Let a be a composition and fi be a partition of the same size. Then, 



pe-P(a) 



(36) 



This conjecture contains Conjecture 5.8 of [I7j and Theorem 17.11 above as special 
cases r = 1 and r = oo, respectively. Also, extensions for paths of more general 
representations without partition /i are discussed in Section 5.3 of [T7] . 

Example 7.3 Let us consider case a = (4, 1, 1) and /i = (4,2). The following is a 
list of paths p and the corresponding value of tau function 2) (p)- For example, the 
top left corner 1 111123 l| means p = |T]®|T]®|T|(8)|T|(8)[2](8)[3]and r^^^ 2) (p) = 1- 



111123 1 


111132 2 


111213 2 


111231 3 


111312 2 


111321 1 


112113 3 


112131 4 


112311 3 


113112 3 


113121 2 


113211 2 


121113 4 


121131 5 


121311 4 


123111 5 


131112 4 


131121 3 


131211 4 


132111 3 


211113 1 


211131 2 


211311 1 


213111 2 


231111 3 


311112 5 


311121 4 


311211 5 


312111 6 


321111 4 



Summing up, LHS of Eg. ( 1361) is 



pGP((4,l,l)) 



which coincides with the RHS of Eq. (l36|) . Compare this with t^_^ 2) data for the same 
set of paths at Example 5.9 of [17]. 



Acknowledgements: The work of RS is supported by Grant-in-Aid for Scientific 
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